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Among the monotone metrics on the (n^ — l)-dimensionaI convex set oi n x n density matrices, 
as Petz and Sudar have recently elaborated, there are a minimal (Bures) and a maximal one. We 
examine the proposition that it is physically meaningful to treat the volume elements of these 
metrics as densities-of-states for thermodynamic purposes. In the n — 2 (spin- 5) case, use of the 
maximal monotone metric, in fact, does lead to the adoption of the Langevin (and noi the Brillouin) 
function — thus, completely conforming with a recent probabilistic argument of Lavenda. Brody and 
Hughston also arrived at the Langevin function in an analysis based on the Fubtni-Study metric. It is 
a matter of some interest, however, that in the first (subsequently modified) version of their paper, 
they had reported a different result, one fully consistent with the alternative use of the minimal 
monotone metric. In this part I of our investigation, we first study scenarios involving partially 
entangled spin-i particles (n = 4, 6, . . .), and then a certain three-level extension of the two-level 
systems. In part II, we examine, in full generality, and with some limited analytical success, the 
cases n = 3 and 4. 
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I. INTRODUCTION 



The Horodecki's have recently studied the relationship between quantum entanglement and the Jaynes inference 
scheme based on the maximization of entropy [^-H . They found that the Jaynes methodology "can lead to problems 
with processing of entanglement" . In this study, we also examine connections between the entanglement of quantum 
systems and their statistical/thermodynamic properties, but from a rather different viewpoint (cf. [^). A principal 
motivation for us in undertaking this work has been the recent contention of Brody and Hughston m that the 
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conventional (Jaynes) density-matrix approach to the canonical ensemble is semiclassical in certain respects (since 
it eliminates the weighting of the quantum phase space volume) , along with their accompanying presentation of an 
alternative "quantum canonical ensemble" based on the metrical geometry of this space. In particular, Brody and 
Hughston analyzed the case of a spin-^ particle in a heat bath and arrived (in the first version ||] of their paper 0|) 
at the partition function, 

QiP) = 2V^r{^)i0h)-'h{(3h) = TT{(3h)-H,{(3h). (1) 

Here, (3 — is the inverse temperature (T being temperature), k is Boltzmann's constant, h is the magnetic field, 
and Ii{x) is a particular modified (hyperbolic) Bessel function {I^{x)^v = 1) of the first kind. (Such functions often 
appear in the distribution of spherical and directional random variables [^.) From (0), one can derive the expectation 
of the energy eq. (14)], 

d\ogQ{(3) ^ hhifih) ^ hif^B/kT) 

dp h{(3h) h{fiB/kT)' ^' 

where ^ is the particle's magnetic moment, and B is the external magnetic field strength, with h = fiB. (Ratios of 
modified Bessel functions, such as appear in (^, play "an important role in Bayesian analysis" In the limit /? — > 0, 
the expected value of the energy (^) is 0, while the variance about the expected value is, then, ■^.) The semiclassical 
analogue of (||) is the Brillouin function ||^, 

Ii (Bh) 

E ^ -h^—— ^ -htanh(3h. (3) 
/_i(/5/i) 

Not only have Brody and Hughston expressed certain reservations and qualifications regarding the Brillouin function, 
so has Lavenda p. 193]. He argued that the "Brillouin function has to coincide with the first moment of the 
distribution [for a two- level system having probabilities and ^^'^JL^ , where x = j^], and this means that the 

generating function is Z{x) = coshx. Now, it will be appreciated that this function can not be written as a definite 
integral, such as 

zm^lJj^^dx^^-^^Qh.m, (4) 

because the integral form for the hyperbolic Bessel function, 

lAx) = . /' e±-*sini-^df, (5) 

V7rr(2 + 1^) J-1 

exists only for v > ^. This means that /_i(a;) = (^)^ cannot be expressed in the above integral form. Since the 
generating function cannot be derived as the Laplace transform of a prior probability density, it casts serious doubts 
on the probabilistic foundations of the Brillouin function. In other words, any putative expression for the generating 
function must be compatible with the underlying probabilistic structure; that is, it must be able to be represented 
as the Laplace transform of a prior probability density" (see also the further detailed remarks of Lavenda, to much 
the same effect |10, pp. 20, 198]). Since the model yielding (|l|) and (j^) is based on the integral form for I^{x), for 



= 1 > i, it clearly accords with the requirements of Lavenda. 

Brody and Hughston ]0] raised the possibility that the quantum canonical ensemble could be distinguished from 
the conventional canonical ensemble by a suitable measurement on a sufhciently small quantum mechanical system. 
In such a case, they argued there would not seem to be any a priori reason for adopting the semiclassical mixed state 
approximation, which allows random phases to be averaged over. Park and Band, in an extended series of papers p]] , 
expressed various qualms regarding the conventional (Jaynes) approach. Park |l^ ] himself later wrote that "the details 
of quantum thermodynamics are presently unknown" and " perhaps there is more to the concept of thermodynamic 
equilibrium than can be captured in the canonical density operator itself." Friedman and Shimony ]ll3|] , in a classical 
rather than quantum context, claim to have "exhibited an anomaly in Jaynes' maximum entropy prescription". 

The (initially presented) results of Brody and Hughston ]|^, that is, (|l|) and, implicitly, (^, had, in fact, been 
reported somewhat earlier by Slater ]fi4| , along with parallel formulas for the quaternionic two-level systems, in 
particular, the partition function (cf. (H)), 
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Q(/3) - 4^/^^(^)(/3/^)-2/2(/3/^) = MPh)-^ HPh) / ^- (6) 

This other analysis [Q, similarly to relied upon a metrical geometry, but the two approaches pursued appear, at 
least superficially, to be somewhat different. The work of Brody and Hughston employed the Fubini-Study metric on 
the complex projective space CP" (the space of rays — which they regarded as the "true 'state space' of quantum 
mechanics"). The study of Slater, on the other hand, utilized the Bures metric, which is defined on the space of 
density matrices ||l7|-[l9|. However, Petz and Sudar have recently shown that the extension of the Bures metric to 
the pure states is exactly the Fubini-Study metric, and that, in point of fact, the Bures metric is the only monotone 
metric which admits such an extension. (The Bures metric is the minimal monotone metric jl5[| . The Bures distance 
— known also as the Hellinger or Kakutani inner product — is defined as an inner product of two probability measures 
on some measurable space [Q. Therefore, it can be naturally defined on the space of pure states, namely, CP".) 
So, these two approaches may be demonstrably fully consistent with one another — as their agreement in yielding 
the results (|l|) and (^) might lead one to hypothesize. However, we must bear in mind that Brody and Hughston 
were led, in the second version of their paper H, to revise certain of their original conclusions H, in the manner 



indicated at the beginning of sec. IV. (In their several recent joint papers dealing with quantum statistical issues, 
Brody and Hughston have chosen to "emphasize the role of the space of pure quantum states, since in the Hilbert 
space based classical-quantum statistical correspondence this is the state space that arises as the immediate object 
of interest. In fact, the space of density matrices has a very complicated structure, owing essentially to the various 
levels of 'degeneracy' a density matrix can possess, and the relation of these levels to one another" [l6|.) 

In the present communication, we follow in sec. || the specific line of argument of Slater (based on the Bures 
metric), with the objective of developing "quantum canonical ensembles" for higher-dimensional situations than 
the (two-dimensional) one presented by a single spin-i particle, previously analyzed. We examine several different 
scenarios for partially entangled spin-i particles, and conduct an analysis of the same nature in sec. HI for a certain 



i/iree-level extension of the two-level systems, previously studied ]2l| ]. Then, in sec. IV, we reexamine the scenarios 
studied in sec. |l|, but now with the use of the maximal rather than minimal monotone metric. In an appendix, we 
report information gains from possible outcomes of joint measurements of the partially entangled spin-^ particles. 

We would also like to bring to the reader's attention, the second part of this study j2^, in which certain results 
(indicated in sec. 0, based on the maximal monotone metric, the subject of sec. IV) of a surprisingly simple nature 



have been obtained for the most general (unrestricted) form of scenario for a spin-1 particle (describable by a 3 x 3 
density matrix with eight free parameters). In ||2^ , we also derived some results pertaining to the (yet unrealized) 
possibility of extending the form of analysis taken there to general spin-|, or equivalently arbitrarily coupled or 
entangled spin-i systems (describable by a 4 x 4 density matrix with fifteen free parameters). (We are confronted 
there with a challenging problem of diagonalizing a 6 x 6 symmetric matrix.) 



II. ANALYSES BASED ON THE MINIMAL MONOTONE (BURES) METRIC 

In general, the 4x4 density matrix (p^^^^^) of a pair of (arbitrarily entangled) spin-i particles (a, 6) can be written 
in the form (we adopt the notation of p^), 

where and cr(°)^(^) are Pauli matrices acting in the space of particle a and 6, respectively. The three-vectors 

^M,ib)^ where = (^i''\ 4"^), correspond (in the case of photons) to the Stokes vectors, while the parameters 
Cy describe the interparticle correlations. If the two particles are independent (nonentangled), then, Qj = ^i'^^^f^- 
In all the scenarios considered below, it is assumed that particle a is described by the same 2x2 density matrix as 
particle b, that is, p'-"-' = p^''\ or equivalently, ^"^"^ = 



A. Particles a and b are Polarized and Correlated with Respect to the Same Direction 

For the first of several scenarios to be examined, let us set twelve of the fifteen parameters in the expansion (|^) ab 
initio to zero, leaving only: (1) which we equate to ^j"''; and (2) Cn. This corresponds to a situation in which 
the two particles are unpolarized and uncorrelated except in one direction (associated with the index "1") of three 
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underlying orthogonal directions. Thus, we are concerned with a rfou6/?/-parameterized set of 4 x 4 density matrices. 
These have the eigenvalues, 

i(l + 2Cj'^)+Cn), i(l-2ei")+Cii) (8) 

and corresponding eigenvectors, 

(0,0,1,0), (0,1,0,0), (1,0,0,0)), (0,0,0,1) (9) 

(Such a system is in a state of degeneracy if either Cii = 1 or Cn = — 1 ± 2^[°'\ since at least one of the eigenvalues 
is, then, zero. For the system to be in a pure state, we must have Cn = 1 and = ±1.) We employed these results 
((H), (H)) in the formula of Hiibner jl^ for the Bures metric, 



dB[P,P + dp) = : (10) 



where is the i-th eigenvalue and < i\ the corresponding eigenvector of an n x n density matrix p. For the case at 
hand, we have the result, 

di3(p,p + dp)2 = 5^(„)^(a)d^J°^ +5^(a)^^^d^J'''dCii +gciiCiidCn, (H) 

where 

^ ' 2(-4e('^) +(1 + Cii)2) 



and 



' Ai^^^ -(1 + Cii)^ 



^ 1 - 2g|°^' + Cll 

4(-i + 2ej'^)-Cii)(-i + Cii)(i + 2ei'^^ + Cii)' 



The corresponding volume element of the Bures metric is (cf. [p4|) 



- '-^^^ - iTiV (-i + 2ei"^-cii)(-i + cii)(i + 2ei") + cii)' ^''^ 

If we integrate this volume element, first, over ^j;"-* from — 5 — to 5 + %" and, then, over Cii from -1 to 1, we obtain 
the result ^. (These limits define the domain of feasible values of and Cii — which determine a triangular region 
— for our doubly-parameterized density matrix. Outside this region, not all the eigenvalues of p^°-'^^ lie between 
and 1, as they must.) Dividing (|l5| ) by ^, we obtain a (prior) probability distribution |2^ ] over the domain of 
these doubly-parameterized 4x4 density matrices. Again, integrating the resultant probability distribution over 
between the same limits as before, we obtain a univariate probability distribution, 

. log(l-fCii)-log(-l-Cii) ^g. 
27r\/2VW^ 

over the interval Cii G [^Ijl]- If now regard (|l^) as a (normalized) structure function or density-of-states for 
thermodynamic purposes, multiply it by a Boltzmann factor, e"'^'"-", and integrate over Cii from -1 to 1, we obtain 
the partition function, 

^ e-^-V^erfi(y2^) 
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where erfi denotes the imaginary error function (The error function erf(z) is the integral of the Gaussian 

distribution, that is, e~*^^dt.) For the expected value of the "energy" (/iCii), we have, then, 

d(3 2/3 erfi(^/2;g7I)^^^• ^ ' 

As /3 — ^ 0, this expected value approaches ^, while the variance about the expected value approaches 

If we choose in this scenario to first integrate the volume element ( |l5| ) over Cn (rather than ^j^"'') from —1 — 2^^"'' 
to — 1 + 2^[°'\ then we obtain (where K represents the complete elliptic integral of the first kind), 

+ (19) 



2V-^i"^-l 



which can not be explicitly integrated over ^^["''^ from -1 to 1 (nor if, first, multiplied by a Boltzmann factor, e '"'^i 

B. Particles a and b are Polarized and Uncorrelated in One Direction, Unpolarized and Correlated in An 

Orthogonal Direction 



Let us modify the previous scenario (sec II A) slightly by now setting the formerly free parameter Cn to 0, and 



letting C22 be free instead. We still maintain = ^[°'\ with the remaining twelve parameters once again set to zero. 
This corresponds to a situation in which the particles a and b are correlated in one particular direction (labeled "2" ) , 
but unpolarized in that direction. The elements of the Bures metric are, then, 

= , (20) 

= , (21) 



and 



_ l-Cl2 + 2g|°^'(-2 + Cl2) ,,,, 

5C22C22 - --2 • l^^i 

4(-l + Cl2)(-l + 4d"' +CI2) 

The integrations of the corresponding volume element of the Bures metric ((pO|)-(p^)) are now performed, first, over 

Ci ' from - V±_i22 v±_^ ^^^^ ^Yicn, over C22 from -1 to 1. (The feasible values lie within an ellipse, the equation 

of which is 4^^"^ + CI2 = !•) This gives us a result of which we can use to obtain a normalized volume element, 
that is, a (prior) probability distribution, 

^ =, (23) 



^Vl-4d"^'-Cl2 



over the domain of the doubly-parameterized 4x4 density matrices for this scenario. The univariate marginal 
distribution of (23) over (22 & ["Ij 1] is simply uniform (|) — which we take as our (normalized) structure function. 



Applying the Boltzmann factor, e f^'^i^^ to it, gives us a partition function (cf. 



yielding an expected value of the "energy" (/1C22) equal to 
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(In the limit /3 0, this expected value approaches 0, while the corresponding variance approaches ^.) It should be 
noted that the results ( p^ and (^) are formally equivalent to those given by the Langevin model of paramagnetism 



C. Particles a and b are Polarized and Uncorrelated in One Direction, Unpolarized and Equally Correlated in 

the Two Orthogonal Directions 



The only difference between this scenario and the previous one (sec. II B ) is that the correlation C33 is not set to 
0, but rather equated to ^22- (So, although the particles a and b are unpolarized in the third direction, the outcomes 
of their individual spin measurements in this direction may be correlated with one another.) This leads to a more 
simple outcome. The normalized form of the Bures volume element is, now, 

\ (26) 

The domain of feasible values is the square defined by the lines ^J"'' = ±| and C22 = (So, the bivariate probability 
distribution ( ^ ) factors into the product of two univariate probability distributions.) We can, then, multiply ( ^6|) by 
the bivariate Boltzmann factor e^^^^^^ '~''c'»C22 g^j^(j integrate over the square region to obtain the (product) partition 
function, 

QiA^Pc) = Ioi^)U^)- (27) 

We have also obtained partition functions of the form /o(^) in two quite different scenarios, in which we set 
thirteen of the fifteen parameters in the expansion (|^) to zero and, otherwise, set the cross-correlation (21 equal to 
either ^12 or to — Ci2- 



D. Particles a and b are Unpolarized, but Equally Correlated in Three Orthogonal Directions 

In this scenario, we set the six components of the two vectors and i^'^''-' to zero, so the particles a and b 
are assumed to be unpolarized in each of the three orthogonal directions (and, thus, with respect to any arbitrary 
orientation). We also fix Qj = if i 7^ j, so there are no correlations between spin measurements in different directions. 
(So, in these respects, a and b are independent or nonentangied.) Finally, we set ^33 — C22 = Ciij so correlations 
are allowed between the measurements of a and b in the same direction. Thus, we are concerned here, not with a 
doubly-parameterized family as in the first two analyses, but with a singly-parameterized one. 

We obtain as our prior probability distribution over the feasible range (.11 S 1/3], 

^ (28) 



7rvi-3CiiyrTciT 

from which one obtains the partition function. 



g(/3)=e^/o(^), (29) 



and an expected value of the "energy" (/iCii) of 



dlogQiP) _h 2hC-f) 

df3 1^^- ^ ' 



So, we again encounter a ratio of modified Bessel functions (cf. (||)) The value of ( |30| ) at /3 = is — while that 

hi 
9 • 



of the associated variance is the square of this, that is ^'^ 
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E. Three and More Unpolarized Particles having Equal Highest-Order Intradirectional Correlations 



The first of several scenarios discussed in this section might be regarded as a i/iree-particlc (a, b, c) analogue of the 
two-particle one presented in the immediately preceding section ( [ID ). For a general 8x8 density matrix (p"''''^) 



representing the joint state of the three particles, we have an expansion analogous to (0). In this expansion, we consider 
all the 63 parameters to equal 0, except for the three (Cm, C222, Csss) representing the highest-order intradirectional 
correlations. We regard these three parameters as having a common value, which is designated Cm- In other words, 
there is a possibly nonzero correlation between the outcomes of spin measurements in some fixed direction for the 
three particles. 

The normalized volume element of the corresponding Bures metric is, then. 



i?(|)Vl2Ci^n-4 



(31) 



where E represents the complete elliptic integral, and the range of feasible values is Cm G [^1/V3, l/\/3]. However, 
no explicit formula for the partition function was found. 

Let us continue this line of analysis to the /oMr-particle case. Now, there are 255 parameters in the expansion 
analogous to (^. We set 252 of them to 0, and equate both of the correlations C2222 and C3333 to Cim, so again 
we are concerned with a singly-parameterized family of density matrices. The feasible range of Ciiii is the interval 
[-1/3,1]. We are able to normalize the volume element of the Bures metric over this interval, obtaining the probability 
distribution (cf. (|8|)), 

^ (32) 



T^V^ - CiiiiVi + 3Ciii 



and the partition function (cf. (29)) 



f3h 2dh 

Q(/3) = e--/o(^), (33) 
giving an expected value of the "energy" (/iCiiii) (cf. (|30|)), 

_ d\ogQ{/3) h 2hC-f) 

dp 3^' loC-f)'' ^ ' 

For /3 = 0, this expected value equals |, while the associated variance is 

For the /iwe-particle analogue, the thermodynamic behavior was precisely the same as for the three-particle case 
discussed just above (pl|), with the replacement of Cm b y Cii iii- For the sia;-particle analogue, the same results 
-(^) were obtained as in the fwo-particle case of sec. II D| (replacing Cii by Ciiim)- 



F. Three Unpolarized Particles having Equal Second Highest-Order Intradirectional Correlations 



We modify the scenarios of the previous section (HE) by, now, requiring the highest-order correlations to be 0, 
while equating all the second-order correlations to each other, obtaining, thereby, the one free parameter. In the 
three particle case, there are nine such correlations — the assumed common value of which, we denote by Ciio- The 
corresponding prior probability distribution over the feasible range Ciio G [^1/3, 1/3] is, then, 

-(4-36Cfio)' ^^^^ 

yielding a partition function. 



QW^Ioi^)- (36) 



The expected value of the "energy" (ftCiio) is, 
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dp 3/o(f)- ^ ^ 

.2 

For /? = 0, this is equal to 0, with a corresponding variance of j^. 

In the four-particle analogue, we have twelve second-highest order correlations — the assumed common value of 
which, denoted Ciiio, has a feasible range of [~ 4^3 1' However, we have been unable to determine the set of 
eigenvectors to employ in the formula ([lO|), and, thereby, can not further pursue the analysis. Interestingly though, 
we have been able to determine the eigenvalues and eigenvectors for the five-particle analogous scenario — in which 
the single parameter Ciiiio must lie in the interval [—\, \]- The volume element of the Bures metric is, then, 

2-\/^l + 3CiiiioVl + 3Ciiiio\/— 1 + 5C11110VI + 7Ciiiio 
We have been unable, however, to either normalize this and/or derive a corresponding partition function. 



III. QUANTUM CANONICAL ENSEMBLE FOR A CERTAIN THREE-LEVEL EXTENSION OF THE 

TWO-LEVEL SYSTEMS 

In our final analysis, we apply the same line of reasoning utilized in the previous scenarios to recent results ]2l| ] 
concerning a particular i/iree-level extension of the too-level systems. These were given by density matrices of the 
form, 



V + z X — iy 

2 - 2w 
X + iy V — z 



P^l\ 2 - 2« I , (39) 



so for V = 1, the middle level is inaccessible and we recover the two-level systems. The normalized volume element of 
the associated Bures metric has been found to be eq. (17)], 

(40) 



Air'^vy/l — V\J v"^ — x"^ — y'^ — 

From this, one can obtain (using spherical coordinates in the integrations), the univariate marginal distribution (an 
asymmetric beta distribution) over the interval v e [0, 1] |^, Fig. 3], 

(41) 



Interpreting this as the appropriate (normalized form of the) structure function, multiplying by the Boltzmann factor 
^-f3hv ^j^j integrating the product over v from to 1, we obtain the corresponding partition function, 

3e-^''((l + 2/jfa)V^erfi(V^) - 2^/^e^'-) 
= ^p^)375 ' (42) 

from which the thermodynamic behavior of an ensemble of such systems (|^ can be deduced. For instance, the 
expected value of the "energy" (hv) is given by 

_ dlogQiP) _ {4l3^h'^ + A(3h + 3)0Ferfi(V^) - 2e^''^/p:(2/3/l + 3) 

dp ~ 2f3i{2(3h + l)0Ferfi(V^) - 2ef>^'^) ' ^ ^ 

As /3 ^ 0, this expected value approaches while the variance about the expected value approaches j^^. 
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IV. ANALYSES BASED ON THE MAXIMAL MONOTONE METRIC 



Subsequent to the distribution of their original preprint D. C. Brody kindly informed me that he and his 
co-author, L. Hughston, had "noticed that the phase space weighting used was incorrect, and as a result we find a 
different expression for the expected energy. Thus, although the claim of the paper is unchanged, we are now revising 
the example considered therein. Interestingly, the same (what we now believe is the correct) result was also predicted 



in our stochastic thermalisation model ( :juant-ph/9711057) — that is, instead of the ratio of the Bessel functions, 



we found the result given by the Langevin function." 

In a second message (also prior to the appearance of 0), Brody wrote: "The only change [in we had to make 
is to replace the 'phase space volume' by 'state density' (i. e. weighted volume). As a consequence, this results 
in changing the calculation involving the spin-^ example. Nevertheless the behaviour of the Langevin function is 
somewhat similar to that of the ratio of the Bessel functions (e. g., energy, heat capacity)." 

This reassessment of Brody and Hughston is paralled by the replacement of the minimal monotone metric by the 
maximal monotone metric (of the left logarithmic derivative), since the use of the latter (coupled with a limiting 
argument), in fact, leads to the Langevin function, that is the ratio and not (For a comprehensive treatment 

of monotone metrics and the "geometries of quantum states", see [p^.) Such a view would, in addition, be appar- 
ently consistent with an earlier study of the author p5[ . There, it was found that the maximal monotone metric 
yielded a more "noninformative" (that is, more "neutral" or less "biased" ) prior, for Bayesian analysis, than did the 
minimal monotone metric (cf. pTf). (In the thermodynamic context, we form posterior distributions by multiplying 
prior distributions by Boltzmann factors, while in the Bayesian context, posteriors are formed by multiplying prior 
distributions by the likelihoods of a certain finite set of outcomes.) 

Let us elaborate somewhat upon these points. The volume elements of both the minimal and maximal monotone 
metrics over the three-dimensional convex set ("Bloch sphere") of spin-^ systems are proportional to expressions of 
the form, 

(44) 



where u = i in the minimal case, and | in the maximal one |T^, eq. (3.17)]. In the minimal case, (^) can be 
directly normalized over the Bloch sphere (Ci +^2 +^1 ^ 1)j ^-^d we are able to obtain univariate marginal probability 



distributions of the type. 



(45) 



Employing this result in Poisson's integral representation of the modified spherical Bessel functions (cylinder functions 



of half integral order) [g^ 29 1 , we recover the associated partition function (|l|) , as originally reported by Brody and 
Hughston 1^, and independently. Slater iQ. For u = |, (^J) is not normalizable (that is, improper) over the Bloch 
sphere. However, it can be normalized over a three-dimensional ball of radius R < 1. Then, integrating out one of 
the three Cartesian coordinates (say, ^3) over this ball of radius R, we obtain the bivariate probability distribution. 



iR'-i)^R^-ei-c. 



27rVl-i?2(-l + + + (^2 _ 1) sin-i R) 

In the limit R I, this converges to the bivariate probability distribution, 

1 



(46) 



2 



(47) 



over the unit disk (^^ -t-^| < 1). Its two univariate marginal probability distributions over ^1 and ^2 are simply uniform 
distributions (^) over the interval [-1,1]. They, then, give rise (again, applying Poisson's integral representation) to 



the Langevin function — as in the revision by Brody and Hughston (cf. sec. [IB). It is also of interest to note 
that if we ab initio set one of the coordinates to zero in the volume element for the minimal metric (|4^), then the 
thermodynamic properties in such a conditional case are identical to those found in the maximal monotone analysis. 
So, we have two independent ways of arriving at the revised (and final) result of Brody and Hughston. 

We have, in fact, conducted several analyses parallel to those reported above in the main body of the paper, but 
based on the maximal, rather than minimal monotone metric. (To obtain the specific forms of the metric, we solved 
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— for the two-particle scenarios — 
opportunity to report these results. 



sets of thirty-two linear simultaneous equations [E8^, eq. (4.26)].) We take this 



For the analysis of sec. [I A, we arrived at a maximal monotone metric simply proportional to the minimal monotone 
metric ((p^-(p^) and, hence, the same normalized density-of-states (16) and partition function ([T^). 
For the analysis of sec. [IB, we obtained a volume element of the maximal monotone metric equal to 



V2VC?i+2ei'^^^-l 



(48) 



Its integral over the domain of feasible values diverges, however, and we have not been able to derive results analogous 
to (||) and (|2|). _ 

Similarly, we 



For the analysis of sec. II C, we arrive at precisely the same results as before, that is, (|26|) and (C 



obtain no differences with the analysis of sec. IIP 



For the three-particle analysis of sec. HE, the normalized volume element of the maximal monotone metric is (cf. 



V3 



so that 



Q(/3)=/o(^). 



(49) 



(50) 



(In the minimal monotone or Bures metric analysis, on the other hand, we were unable to find an explicit formula for 
the partition function.) For the four-particle analysis of sec. [I E , the analogous results based on the maximal monotone 
metric are precisely the same ((^2|) - (p^). The maximal- monotone- metric analysis of the five-particle scenario yielded 
the same outcomes ((^9|)-(p0|)) as the maximal-monotone-metric analysis of the three-particle scenario. 

For the three-particle scenario of sec. [IF, an analysis based on the maximal monotone metric yielded the same 
results ((|35|)-(p7|)). As before, we have not been able to compute the volume element for the associated four-particle 
scenario. 



For the analysis of sec. HI, corresponding to a certain three- level extension of the two- level systems, we obtain a 
volume element (cf. (HOh), 



1 



Vl — v{v'^ — — — 2^) ' 



(51) 



not proportional to ([lO]). However, using spherical coordinates and a limiting argument (that is, integrating over the 
radial coordinate r from to normalizing this result, and then letting i? — > 1), it can be shown to yield precisely 
the same univariate marginal probability distribution ( |4l| ) as for the analysis based on the Bures metric, and hence 
lead to the same partition function (^2|). 

It would be of interest to attempt to find a characterization of those cases in which the maximal and monotone 
metrics are the same (and so, of course, are their volume elements) or, more generally (as in the three-level extension 
example), of the cases in which the metrics may not be identical, but the associated thermodynamic properties are, 
nevertheless, the same. 

Let us further note that for the /iwe-dimensional convex set of quaternionic two-level quantum systems, the volume 
elements of the maximal and minimal monotone metrics are proportional to expressions of the form (cf. (^), ]3C| , 
eqs. (7),(20)]), 



(52) 



^ in the minimal case and § in the maximal case. We then find — again applying Poisson's integral 



where u — ^ u^^^ ^^^^^^^^^^^^ ^^^^ ^^^^ 2 
representation to the univariate marginal probability distributions obtained from ( |52D — that in the minimal case, 
the expected energy is proportional to the ratio ^ and in the maximal case (applying a similar limiting argument, as 



used to obtain (47)) to the ratio -p-^. So, these results are quite analogous to those for the three-dimensional convex 



set of complex two-level quantum systems. 
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Let us note that the results for the minimal and maximal monotone metrics coincide for precisely those scenarios in 
which we are concerned with a mutually commuting set of density matrices. Then, the scenario is essentially classical, 
rather than quantum, in nature. Hence, the monotone metric — the one associated with the Fisher information — 
is simply unique. We will also point out that in a related analysis p2[ , we found the maximal monotone metric to be 
considerably more amenable to our (MATHEMATICA) computations than was the minimal monotone metric. 



V. CONCLUDING REMARKS 

As with the instance of a single spin-i particle studied in — giving the previously reported results (Q) and (H) 

— the possible applicability (to small systems, in particular ra]) of the several quantum canonical ensembles presented 
in this letter, awaits experimental examination. (We should note that it is usually considered that entangled particles 
will decohere before they are thermalized, that is, entanglement will be lost on time scales short compared to those of 



thermal relaxation processes associated with energy exchange with the bath |34|.) 

Although wc have examined a number of possible scenarios here, there is clearly much opportunity for further 
systematic exploration along related lines. In fact, in |Q, we were able to assign a prior probability distribution of 
the form ^^^^^^^^ to the two-dimensional simplex of diagonal entries — a,b,c — of the eight-dimensional convex 
set of spin-1 density matrices. (The evident symmetry between these three entries had been broken by a particular 
transformation, suggested by work of Bloore |^^.) Then, we obtained a spin-1 counterpart — Fig. 4] — to the 
Langevin function (7^). These results pertained not to the minimal monotone metric (which proved to be more 



computationally problematical), but only to the maximal one — which is the subject of sec. IV — and is consistent 
with the reassessment by Brody and Hughston of their original work |^ , as the maximal monotone metric applied 
to a single spin-^ system, yields the Langevin function. As noted above, this also accords with certain arguments 
of Lavenda ||l^, pp. 20, 193, 198], who contends that the Langevin function has a sound probabilistic basis, while 
the Brillouin function does not. Also, as previously pointed out, the maximal monotone metric has been shown in 
two separate analyses p5|j27| to be more noninformative in character than the minimal monotone metric. Other 
information-theoretic properties of these metrics have been studied in (cf. ||3^]). 

This study might be viewed as an attempt to respond to the concluding admonition of Petz and Sudar ||l^, p. 2672] 
that "more than one privileged metric shows up in quantum mechanics. The exact clarification of this point requires 
and is worth further studies" . In particular, we have been concerned with evaluating the various thermodynamic 
implications of these monotone metrics. 
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APPENDIX. Information Gains from Differing Possible Outcomes of Joint Spin Measurements 

We avail ourselves of the several prior probability distributions (normalized volume elements of Bures metrics) 
reported in sec. ^ to obtain the information gains for various outcomes of joint measurements of the partiaUy 
entangled spin-i particles, (we might also pursue similar analyses using the results of sec. IV, based on the maximal 



monotone metric.) We pass, in this manner, from thermodynamic considerations to ones of a fundamentally Bayesian 
nature. We follow a methodology previously employed for unentangled spin-^ particles ||35|| . 

We associate with each specific outcome of a joint measurement of the partially entangled particles, a likelihood 
function, the product of which with the prior distribution can be normalized to yield — via Bayes' rule — a posterior 
probability distribution. We then compute the information gain (KuUback-Leibler statistic) of the posterior with 
respect to the prior. The products of the gains and their corresponding likelihoods averaged over all the possible 
outcomes with respect to the prior, then, give the expected information gain. 



For the scenario of sec. II A, if we measure the spins of particles a and h along the axis of polarization, and find them 
to disagree (that is, one spin "up" and the other, "down"), we gain log 3 — | ~ .431946 nats of information, while if 
they agree (both spins either up or down), we gain considerably less, that is, log 6 — f ~ .125093. (The likelihood to 

be used in the Bayes' rule of obtaining a disagreement is ^"^"^^^^ , while that for an agreement is ^"^^^"^ . The expected 
information gain is .329662.) If after one agreement, we obtain a disagreement, then the information gain due to the 
second measurement is greater than that from an initial disagreement, that is, log 5 — y| « .542771. If, instead, there 
is a second agreement, the additional information gain due to that outcome is only .0427712 nats. Further, if an initial 
disagreement is followed by another, the gain is log | — | « .110826, while if the initial disagreement is succeeded by 
an agreement, the gain is log 10 — « .235918. 
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If we were to alter the scenario of sec. II A , setting S^^^ equal to the negative of ^j"' rather than to (}^^ itself (so 
that the spins would be anticorrelated), we anticipate obtaining the same expected information gains, but for the 
reversed situations. For example, the gain due to an initial agreement would be .431946 nats, and not .125093, which 
would be that associated with an initial disagreement, .... 
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Contrastingly, for the scenario of sec. IIB, the information gain from an initial agreement is the same as that for 
an initial disagreement, that is log 2 — ^ ~ .193147, so this is also the expected gain. If a second measurement yields 
an outcome different from the first, the information gain from this measurement is log 3 — | ~ .265279, while if the 
outcome of the second measurement is the same as the first, the gain is log I - i « .0721318 nats. 

For the scenario of sec. [ID, an initial agreement of spin measurements of particles a and b along an arbitrary axis 
yields an information gain of .306853 nats (saturating the Holevo-type bound ||3|,^), while an initial disagreement 
gives .0646381 nats. (The expected information gain is .145376 nats.) A second agreement yields .0680544 nats, while 
a second disagreement gives .0470689 nats. An agreement after an initial disagreement results in a gain of .427868 
nats, while a disagreement after an initial agreement provides .0516789 nats. 

There, of course, still remains the possibility of using positive-operator- valued measures (POVM's) and/or stan- 
dard measurements on multiple copies of (partially) entangled systems of spin-i particles ||35|-|37[|. Along sim- 
ilar lines, in |3^ it was found (based on the corresponding Bures metric) that for a pair of identically pre- 
pared unentangled spin-i particles, one could expect to gain: (a) .357229 nats of information with the use 
of a particular POVM proposed by Peres (having a continuum of possible outcomes); (b) .313478 nats 
for a certain joint measurement scheme of Popescu and Massar ||3^] (cf. |^,Q), having four possible out- 
comes; (c) .280372 nats if the spins of the two particles were to be measured separately in orthogonal di- 
rections; and (d) .212642 nats if measured separately in the same direction. A spin measurement on a 
single spin-i particle can be expected to yield .140186 = •'^^^^^^ nats M. 
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